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Abstract 

The method of non-perturbative renormalization group (NPRG) is apphed to the analysis of 

dynamical chiral symmetry breaking (D^SB) in QCD. We show that the D^SB solution of the 
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^SJ ■ NPRG flow equation can be obtained without the bosonization. The solution, having the singular 

$—1 . 

pi |. point, can be authorized as the weak solution of partial differential equation, and can be easily 
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evaluated using the method of the characteristic curve. Also we show that our non-ladder extended 
approximation improves almost perfectly the gauge dependence of the chiral condensates. 
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I. INTRODUCTION 

The approach of the "non-perturbative renormahzation group" (NPRG) based on Wil- 
son's idea allows us to analyze dynamical chiral symmetry breaking (D^SB) in QCD with 
systematic approximations compared to the ladder Schwinger-Dyson analysis which suffers 
from the strong gauge dependence. The NPRG approach has been formulated as several 
forms of the functional differential equation, so people often call these formulations the 
"functional renormahzation group". Among them, we adopt the flow equation derived by 
Wetterich (ij. 

In NPRG, the 4-fermi coupling constant dynamically generated by the gauge interaction 
blows up at an infrared scale as a signal of D^SB. Due to its explosive behavior, the access 
to the broken symmetry scale below the critical scale is a nontrivial problem. In the ear- 
lier studies, the auxiliary field (bosonization) |3|, |j] , the dynamical bosonization[5j], and the 
explicit breaking term such as the bare mass term |6| have been introduced to access the 
broken symmetry scale. However, using the weak solution method^, we can directly access 
the broken symmetry scale only with the multi-fermi interactions but without any help by 
those methods. 

II. TRUNCATED EFFECTIVE ACTION 

In order to study the dynamical chiral symmetry breaking in QCD with three-flavor 
massless quarks, we define the following truncated effective action: 



^Am=l[^F;:^Fr + ^id,A^)' + ^iZ^^ + tgs4)i^-Vi^,^;A)^, 



(1) 



where we use the covariant gauge with the gauge- fixing parameter ^, and do not represent 
the ghost sector for simplicity. In this definition, the multi-fermi interactions are represented 
by V{tp,il!]A), which we call the fermion potential. 

We substitute the effective action ([1]) into the Wetterich flow equation, and expand it in 
terms of the fields and the derivatives, so that we can obtain the coupled non-perturbative 
RG equations of these coupling constants. Then, lowering the cutoff scale A by solving the 



^ The authors greatly appreciate helpful comments by Prof. Akitaka Matsumura who told us how to 
construct the weak solution. 



RG equations, the gauge interaction dynamically induces multi-fernii interactions. Among 
them, the scalar 4-fermi interaction brings about the D^SB at an intermediate cutoff scale 
as the Nambu-Jona-Lasinio model does. The 4-fermi coupling constant blows up at a finite 
infrared scale Ac, showing the signal of the D^SB, and consequently its RG flow cannot go 
beyond the critical scale Ac. 

III. PARTIAL DIFFERENTIAL EQUATION 

To go beyond the critical scale Ac, we solve the flow equation without expanding the 
fermion potential with respect to the scalar fermion-bilinear field a{= ipip) that is the central 
operator for D^SB. For simplicity, the arguments of the fermion potential are limited to the 
scalar operator: V{tp,il!]A) — t- V{a]A). Besides, the interactions are limited to the ladder 
type ones, and consequently we obtain the ladder flow equation as the following partial 
differential equation (PDE) J3|, \j^ : 



A* 
dtV{a; t) = -rj^ad^V + —\n 



,2 N 2 



, C2 glcy 
4A2 



(2) 



where C2 is the second Casimir invariant of S'f/(3). Here the dimensionless scale parameter 
t{= logAo/A) and the anomalous dimension of the quark field rj^{= dt\ogZ^) are intro- 
duced. The PDE 02]) can be numerically solved by replacing the derivatives with the finite 
differences. Practically, we work with the PDE for the first derivative of the fermion poten- 
tial, M(o";t)(= daV), which we call the mass function. As shown in Fig. [H the numerical 
solution of the mass function at cr = has a finite jump below the critical scale Ac. This 
is nothing but the signal of D^SB, and the infrared value of the mass function at the limit, 
(J —7- +0, is nothing but the dynamical mass mdyn. 

IV. WEAK SOLUTION 

While the singular behavior have numerically been evaluated, the non-analytic solution 
cannot be mathematically authorized as a solution of the PDE. Such a solution can be 
defined as the "weak solution" which satisfies the integral form of the PDE[7|. Here we skip 
the detail of its definition, but briefly introduce the method of the characteristic curve to 
evaluate the weak solution. 
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FIG. 1: RG evolution of the mass function M{a;t). 



Representing the right-hand side of Eq. (j2j) as a function F{a,M;t), the PDE for the 
mass function can be written by 



dtM{a;t) 



-F„ - d^M ■ F. 



Mi 



(3) 



where we use the notations, F„ = d^jF^cr, M; t) | Af , Fm = dpiF^a, M; t)\a-. Now let us consider 
the mass function on a characteristic curve, a = a{t; (Tq), where a^ is the initial position of 
a. If a{t; o"o) obeys the following differential equation. 



d _. , ^ 
— a(t; ao) = Fm, 



(4) 



then the derivative d^M vanishes in the PDE: 



d 



-M(.;i) = -F, 



(5) 



Eqs. dl]) and (|5]) are coupled ordinary differential equations equivalent to the PDE (|3]). 
The parametric curve of the mass function M{a; t) in terms of ctq can be constructed by 
solving them for each ctq as shown in Fig. [2j Below the critical scale Ac, the parametric 
curve is no longer regarded as the one-to-one mass function of a. However, a condition 
of the weak solution called by "Rankine-Hugoniot condition" uniquely gives the one-to-one 
function shown in Fig. [TJ8|. 
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FIG. 2: Parametric curve of the mass function M(a;t) in terms of ctq. 
V. RESULT AND DISCUSSION 

Finally we show the result of the non-ladder extended flow equation, where the non- 
ladder corrections are added to the ladder flow equation (I2])|2[. In this approximation, 
the commutator of the generator of SU{3)c is ignored consistently for obtaining the gauge 
independent result |9|. 

The gauge dependence of the chiral condensates {ipi'), which is actually obtained by 
introducing the bare mass term of quarks as its source term, is shown in Fig. [31 The 
gauge dependence of the non-ladder extended results is greatly suppressed compared with 
the ladder results. Therefore we conclude that our non-ladder extended flow equation well 
respects the gauge independence. Here it should be noted that, in the Landau gauge (.^ = 0), 
the gauge-dependent ladder result of the chiral condensates coincides with the almost gauge- 
independent non-ladder extended one. This feature of the Landau gauge proves a folklore 
that the ladder approximation looks good particularly in the Landau gauge. 
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FIG. 3: The dependence of the chiral condensates (V'V') ^^ the gauge-fixing parameter ^. 
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